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ABSTRACT 

Context. Liberal evolution of interacting binaries has been proposed previously by several authors in order to meet various observed 
binary characteristics better than conservative evolution does. Since Algols are eclipsing binaries the distribution of their orbital 
periods is precisely known. The distribution of their mass ratios contains however more uncertainties. We try to reproduce these two 
distributions theoretically using a liberal scenario in which the gainer star can lose mass into interstellar space as a consequence of its 
rapid rotation and the energy of a hot spot. 

Aims. In a recent paper (Van Rensbergen et al. 2010, A&A) we calculated the liberal evolution of binaries with a B-type primary at 
birth where mass transfer starts during core hydrogen burning of the donor. In this paper we include the cases where mass transfer 
starts during hydrogen shell burning and it is our aim to reproduce the observed distributions of the system parameters of Algol-type 
semi-detached systems. 

Methods. Our calculations reveal the amount of time that an Algol binary lives with a well defined value of mass ratio and orbital 
period. We use these data to simulate the distribution of mass ratios and orbital periods of Algols. 

Results. Binaries with a late B-type initial primary hardly lose any mass whereas those with an early B primary evolve in a non- 
conservative way. Conservative binary evolution predicts only ~ 12 % of Algols with a mass ratio q above 0.4. This value is raised up 
to ~ 17 % using our scenario of liberal evolution, which is still far below the ~ 45 % that is observed. 

Conclusions. Observed orbital periods of Algol binaries larger than one day are faithfully reproduced by our liberal scenario. Mass 
ratios are reproduced better than with conservative evolution, but the resemblance is still poor. 
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1. Introduction 

Warner (1978) introduced the denomination liberal to distin- 
guish binary evolution with a mass loss and subsequent angu- 
lar momentum loss from the conservative case where no mass 
leaves the system. Van Rensbergen et al. (2008) developed a 
liberal scenario in which mass can be lost from a binary dur- 
ing a short era of rapid mass transfer soon after the onset of 
Roche lobe overflow (RLOF). The joint catalog with evolution- 
ary tracks was extended with the evolutionary tracks for binaries 
starting from sufficiently small orbital periods, so that RLOF 
starts during core hydrogen burning of the donor (case A) as 
presented by Van Rensbergen et al. {2010). In this contribution 
these calculations have been extended towards larger initial or- 
bital periods so that RLOF starts after exhaustion of hydrogen in 
the core of the donor (case B). Our catalog (Van Rensbergen et 
al. [2008 ) contains now 561 conservative and liberal evolution- 
ary tracks and is available at the Centre de Donnees Stellaires 
(CDS). The grid of calculations covers binaries with a B-type 
primary at birth. 

A binary is considered to be an Algol when the semi- 
detached system shows the typical characteristics as mentioned 
by Peters (2001), where the less massive donor fills its Roche 
lobe, the more massive gainer does not fill its Roche lobe and 
is still on the main sequence and the donor is the cooler, fainter 
and larger star. The distribution of orbital periods of Algols is 
precisely known since the time between two subsequent eclipses 
for every system has been measured accurately. The catalog of 



Budding et al. (2004]! extended with the semi-detached Algols 
from Brancewicz et al. (1980) supplies us with 303 Algols which 
can be issued from binary evolution with a B-type primary at 
birth. Van Rensbergen et al. (2006) showed that the distribution 
of mass ratios strongly depends on the method used to interpret 
the observations. Substantial differences are found between mass 
ratios derived so as to make the parameters of the most massive 
star fit main sequence characteristics (qMs) and mass ratios qic 
and qs d which are respectively obtained by the light curve so- 
lution and using the fact that an Algol is a semi-detached sys- 
tem. In this paper we used the observed mass ratio distribution 
of Van Rensbergen et al. (2008 ) who applied a mix of qMs and 
qic which represents the numbers in the Pourbaix et al. catalog 
(I1ES best. 

Mass loss during liberal evolution is defined by a quantity f3 
G [0,1] giving the fraction of mass lost by the donor (subscript 
d) that is accreted by the gainer (subscript g): 

M g = -(3Mf 0F (1) 

Conservative evolution is characterized by [3 — 1 . Liberal the- 
oretical calculations (i.e. those in which, at least for some time, 
/3 < 1) depend not only on the amount of mass lost from the 
system (related to the parameter j3) but also on the amount of 
angular momentum taken away by this matter (related to the pa- 
rameter a), formulated by the a, jS-mechanism of Podsiadlowski 
et al. (11992b . The quantity a is then defined by Rappaport et al. 
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(11983b so that, for the angular momentum loss rate from the sys- 
tem / in the case of mass loss rate of the donor M R , L0F one has 

d 

that: 

/ = a Mf^(l-/J)^ (2) 

with a and P respectively the semi-major axis and orbital 
period of the binary system. The value of a thus depends on 
how much angular momentum is carried away from the system 
by a certain amount of mass loss, i.e. on the way in which matter 
leaves the binary. 

Values of B < 1 were proposed by Meurs & Van den Heuvel 
( 1989 ) in order to theoretically reproduce the observed numbers 
of persistent strong massive binary X-ray sources and the num- 
ber of observed Wolf-Rayet binaries. Binary evolutionary calcu- 
lations with constant values of B (e.g. B = 0.5) were then repro- 
duced by various authors, e.g. De Loore & De Greve ( 1992J and 
more recently by De Mink et al. d2007l > (B = 0.25, 0.5 and 0.75). 

Although mass loss from binaries is needed in evolutionary 
theory, it would be astonishing that B does not depend on the 
mass transfer rate. It is more plausible that evolution of semi- 
detached binaries remains conservative (J3 = 1) during the long 
lasting quiet eras of slow RLOF, and that mass can only be lost 
from the system [B < 1) during short lasting violent eras of rapid 
mass transfer. 

The stellar radiation and the spin-up of the gainer as given by 
Packet dl981b is used by Langer ( 1998) to compute the mass loss 
rate for massive binaries. This scenario yields a time dependent 
behavior of B as published by Wellstein et al. (2001 ) and Petrovic 
et al. ( 120051 

The stellar radiation of a rapidly rotating intermediate mass 
gainer is however not sufficient to drive mass out of the system. 
Hence we introduced in the Brussels binary evolution code the 
possibility of mass loss from a binary as a consequence of spin- 
up and hot spots created on the gainer by mass transferred from 
the donor. The physics of this hot spot model, as well as the 
resulting mass loss from a binary with a B-type primary at birth 
were discussed by Van Rensbergen et al. (2008 , 2010) using this 
scenario. 

2. Method 

2. 1 . Brussels evolutionary code 

Our results were obtained using the Brussels binary evolution- 
ary code. A detailed description of this code can be found in De 
Loore & Doom ( I1992| l. It is in fact a solution of the stellar struc- 
ture equations, introduced by Henyey et al. (I19591 1 1964b and 
adapted by Kippenhahn & Weigert (1967) and Kippenhahn et 
al. ( 119671 1 for binary evolution. The calculations were performed 
with the OPAL-opacities (Rogers & Iglesias [l992l Iglesias et al. 
[19921 Iglesias & Rogers [T996I 

For the input physics, convective mixing, radius corrections 
and nuclear physics, we refer to Prantzos et al. 0986). Moderate 
convective core overshooting is applied, as described in De 
Loore & De Greve ( 1992). Mass loss by stellar wind is included 
using the semi-empirical formalism of De Jager et al. 0988), 
and the period changes due to subsequent angular momentum 
loss are also taken into account. However, these small losses are 
of little importance to the evolution outcome. For the spin-up 
of the gainer star and the evolution of the orbital separation, an 
angular momentum balance is considered, which includes the or- 
bital angular momentum of the system and the spin angular mo- 



menteum of the gainer. The latter is calculated assuming solid- 
body rotation. When the gainer is spun up, its increase in spin 
angular momentum is extracted from the orbit. The spin angular 
momentum of the donor, which is small compared to the others, 
is currently not taken into account, but will be in future research. 

2.2. Initial conditions for binaries with a B-type primary 

Van Rensbergen et al. (|2006) used non-evolved systems in the 
9th catalog of Spectroscopic Binaries of Pourbaix et al. (2004) 
to establish the initial conditions for the evolution of binaries 
with a B-type primary at birth. We distinguished between late 
B-type primaries in the mass-range [2.5,7] M Q and early B-type 
primaries in the range [7,16.7] M Q . The following initial condi- 
tions were adopted: 

• future donors follow a normalized IMF: ^(Md)-C M~, x with 
x-235 as given by Salpeter d 1955b . 

• binaries with a B-type primary at birth have a normalized 
bimodal period distribution n(P)=p (with A such that the 
distribution is normalized over the entire period range) for their 
initial orbital periods (Van Rensbergen et al. 2006): 79.3 % of 
the systems with a late B primary follow a distribution as 
proposed by Popova et al. {1982) from P„„„ = 0.93 d to P max = 
33.37 d, whereas the remaining 20.7 % follow a similar 
distribution between P m ,„ = 33.37 d and P max = 9000 d. Out of 
the systems with an early B primary at birth, 79. 1 % follow a 
Popova distribution from P min = 1.02 d to P max = 12.91 d, 
whereas the remaining 20.9 % follow a similar distribution 
between P m - m = 12.91 d and P max = 4000 d. Such a bimodal 
distribution has also been found by Sana et al. (2009) for the 
initial distribution of orbital periods of massive binaries in 
NGC 66 1 1 , where they found 60 % between P m ,„ = 2 d and 
Pmax — 10 d, whereas the remaining 40 % have initial orbital 
periods between P mi „ = 10 d and P max = 3200 d. 

• initial mass ratios (q=^f) follow a normalized distribution: 

y ¥(q)=C (1 + qY 6 ; (5=0.65 for non-evolved binaries with a late 
B primary and 5=1.65 for non-evolved binaries with an early B 
primary, as given by Van Rensbergen et al. (2006). 

2.3. Liberal evolution of binaries 

We completely followed the scenario as outlined in Van 
Rensbergen et al. ( 120101 1. The details of the scenario can be 
found in that paper. During mass transfer the gainer is spun up 
and a hot spot is created in the accretion region which is located 
at the surface of the gainer in the case of direct impact or at the 
edge of the accretion disk if the criterion of Lubow & Shu (119751 1 
tells us that such a disk is created around the gainer. 

2.3.1 . Spin-up of the gainer 

Disregarding tidal interactions, Packet d 1 98 1 b showed that the 
spin-up of the gainer goes very fast during the short era of rapid 
mass transfer soon after the onset of RLOF. Figure [T] shows 
the spin-up for the gainer of a (6+3.6) M G binary with an ini- 
tial orbital period of 3 days. The spin-up was calculated with 
both strong and weak tidal interaction. The formalism was taken 
from Zahn d 1 977b . Strong tidal interaction is calculated with 
fsync =0.1, whereas f sync = 1 implies weak tides, as adopted by 
Langer et al. (f2003) and discussed in detail by Wellstein ( 12001b . 
One sees in Fig. [T] that the gainer is spun up to critical velocity 
very quickly during the rapid era of mass transfer during core 
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hydrogen burning of the donor. This is the case for slow tidal in- 
teraction as well as for a binary undergoing strong tides. During 
the slow phase of mass transfer the rotation of the gainer syn- 
chronizes rapidly in the case of strong tidal interaction. In the 
case of weak tidal interaction, synchronization is never achieved 
before the onset of RLOF B. The rotation of the gainer spins up 
again towards critical velocity when a second era of fast mass 
transfer occurs during hydrogen shell burning of the donor. It 
is clear that with rapid rotation a favorable situation is created 
for the gainer to lose mass into interstellar space. The example 
shown in Fig. [T] will turn out to evolve marginally liberal: no 
mass will be lost from the system during core hydrogen burning 
of the donor whereas only AM = 0.006 M Q will be lost during 
hydrogen shell burning of the donor due to the combined action 
of rapid rotation and the creation of a hot spot on the gainer's 
surface. 



is often used to calculate the 



2.3.2. The hot spot 

The expression L^ cc = G 

accretion luminosity. This expression is only valid in the case of 
matter coming from infinity impacting on the gainer's surface. In 
the case of interacting binaries the transferred matter is however 
issued from the first Lagrangian point L\ and not from infinity. 
Hence, the real accretion luminosity has to be multiplied by a 
factor D (D e [0, 1]): L acc = D L™ cc . The geometrical factor 
D will be large when L\ is far away from the point of impact 
P but can also be very small for binaries which are almost in 
contact. The factor D is the geometrical factor included in the 
last column of TableQ] The factor S is the fractional surface area 
of the hot spot. This factor has been observed in some cases and 
Pringle ( 1985 ) derived a theoretical expressior^jstarting from the 
characteristics (surface and opening angle) of the opening in the 
first Lagrangian point L\ and the fact that the stream keeps its 
shape until it arrives at the point of impact P: 



0.4556 10~ 4 T - 5 



eff,d p 



\LiP\ (R^ 



Md_ 



(3) 



The efficiency of the liberal scenario is further defined by 
the value of the factor of radiative efficiency K which defines the 
fraction of L acc that is used by the hot spot as radiation pressure. 
This accretion luminosity is weakened on the one hand by the 
fact that only a fraction can be converted into radiation (repre- 
sented by a weakening factor A > 1 : L„dd = ^) and on the other 
hand strengthened because the energy of the impacting material 
is concentrated in a hot spot which is significantly smaller than 
the entire gainer's surface (represented by the fractional surface 
(S = < 1). The radiation pressure of the hot spot is thus 
supplied by 



AS 



(4) 



Van Rensbergen et al. (2010 ) showed that in the case of direct 
impact on the gainer's equator, the quantity K can be calculated 
as follows: 



K = 



4 

spot 



eff,g) 



(5770) 4 



(5) 



1 containing as a printing error, which has to be corrected when ap- 
plying this formula, an exponent (-0.5) instead of (0.5) for T e ff. 



There are only 13 reliable hot spot temperatures available in 
the literature. According to the criterion of Lubow & Shu ( 119751 
and as shown in Fig. [2] ten systems are direct impact systems, 
while the three others have a transient accretion disk. A perma- 
nent accretion disk, emitting continuously in Ha, can only be 
formed when the size of the gainer is sufficiently small relative 
to the size of the orbit. 

In the case of the formation of a hot spot on the edge of 
an accretion disk we have to replace in Eq. |5]) R g by Raisk- We 
further replaced T e ff g by Tdisk for /? Lyr which has an opaque 
and optically thick accretion disk. The adopted value of Tdisk ~ 
8900 K will be discussed in Sect. |A.2| Two other systems have a 
transparent and optically thin accretion disk so that their hot spot 
is seen upon the effective temperature of the gainer star. Table 
[l] contains all the data to find a value for K for 13 interacting 
binaries, with a (rather poor) best fit as a function of the total 
mass of the system, which is luckily only slightly different from 
the expression found by Van Rensbergen et al. (1201 Oi l at the time 
when data were available on only 1 1 interacting binaries: 



* = 3.9188 |^ + ^ 



1 .645 



(6) 



This result allows us to calculate the radiative efficiency r\ of 
accretion onto a spinning main sequence gainer with a hot spot. 
Using L a dd = L^ cc D K S — rj M c 2 with the mass-radius relation 
determined from the stellar evolutionary models of Schaller et 
al. ([1992 ) for gainers in the mass range considered in this paper 
and in the middle of their main sequence life, one obtains: 



n = 5.85 10" 



Mo 



0.446 



M d M g 

M e M Q 



DS. 



(7) 



As expected, we find that this radiative efficiency is low as 
compared to the values found for accretion onto white dwarfs (77 
~ 10~ 4 ), neutron stars (77 x 0.2) and black holes (77 up to 0.423 
for the most rapidly rotating black hole). 

2.3.3. Loss of mass 

The effects of spin-up and accretion luminosity induced photon 
pressure increase with higher mass accretion rates. The critical 
mass transfer rate M^ L0F ' cr " is that at which both effects will 
just compensate the local gravity, so that the local Eddington 
luminosity is reached. When the mass loss rate of the donor ex- 
ceeds this critical value, the gainer will nevertheless accept only 

' RLOF c fit 

a rate M d ' , and the surplus will be lost from the system, 
thus resulting in < 1. The exact formulation of ^ L0F ' CI " as 
a function of evolutionary parameters, as included in the code 
and thus used to calculate the value of /? at each moment of the 
evolution, is derived as equations (21) (for direct impact) and 
(22) (for systems with an accretion disk) in Van Rensbergen et 
al. (f2008). These equations imply that each amount of energy 
extracted from the accretion stream can only be used once, i.e. 
either for spinning up the gainer or for developing or maintaining 
its hot spot. It also shows that when f3 < 1, this does not necessar- 
ily mean that critical rotation is achieved, nor that the accretion 
luminosity alone exceeds the Eddington luminosity, but that both 
effects combined exceed the local gravitation. 



2.3.4. Loss of angular momentum 

Our liberal code thus calculates a time dependent /3(t) self- 
consistently within the model as described in the previous sub- 
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Fig. 1. Equatorial velocity of the gainer of a (6+3.6)M Q binary with an initial orbital period of 3 days. Strong tidal interaction is 
represented by a dashed line, weak tides by a long-dashed one. The critical velocity is indicated by a solid line. 
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Fig. 2. (q - r) diagram for the 13 interacting binaries that were used in our analysis. Algols (q < 1) are at the left of the vertical line 
q = 1.0. The ten direct impact systems are located above the curve ojj (filled squares). The three systems with a transient accretion 
disk around the gainer are located between the curves w m ,„ and Ud (filled circles). Permanent accretion disks would be located below 
oJmin- Systems with q > 1 have not reached Algol characteristics yet. 



section, and assumes that matter is lost from the hot spot on the 
gainer (or accretion disk around the gainer) so that the escaping 
matter removes only the angular momentum of the gainer's or- 
bit. However, the near-critical rotational velocity of the gainer 
considerably enhances its spin angular momentum. The influ- 
ence of the latter in the angular momentum loss process has not 
been considered in this study, but will be included in future re- 
search. One thus obtains a time dependent value of a, describing 
the angular momentum loss from the system: 



M d 



K M g + M d 



(8) 



A typical value of a during the fast and liberal era of RLOF 
A is then a w 0.25 (q ~ 1). When this liberal era is succeeded by 
a fast and liberal era of RLOF B the value of a turns out to be 
much smaller and can easily be calculated from relation dHl. 



2.3.5. Values and uncertainties for TableQ] 

Table [TJ does not contain the values of the present-day orbital 
periods of the systems since these values are very well known 
for the eclipsing binaries used in our analysis. The other values 
used to construct Table [T] are discussed in Appendix A. 

Relation (|6]l was used in our code to gauge the radiative ef- 
ficiency of accretion by a main sequence gainer. Figure|3]shows 
the large uncertainty that remains in the determination of K 
caused by the errors and uncertainties on the values quoted in 
Table [T] The curves above and below respectively add and sub- 
tract one standard deviation from the values of K as obtained 
from relation ( 6]). It is clear that binary evolution would be more 
conservative than argued in this paper if the lower curve were 
valid and less conservative for the higher curve. 

Figure [3] shows error bars obtained from the data quoted in 
Appendix A. The calculation of K requires firstly the knowl- 
edge of the orbital period and of the masses, radii and effective 
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Table 1. Data used to determine K for 13 semi-detached binaries. 



System 


M d + M g 


Rci + R g 


T e ff,d 


Teff,g 


dM/dt 


T spot 


K 


D 




(M s ) 


(Rq) 


(K) 


(K) 


(M a yr- [ ) 


(K) 






VW Cep 


0.90 + 0.25 


0.93 + 0.50 


5000 


5200 


8.769E-08 


7076 


6.528 


0.0445 


ARBoo 


0.35 + 0.90 


0.65 + 1.00 


5398 


5100 


1.484E-07 


5539 


1.064 


0.0535 


CN And 


1.30 + 0.51 


1.43 + 0.92 


6500 


5911 


1 o 1 rr 1 m 

1.215E-07 


6485 


8.251 


0.0261 


KZ> rav 


n on i i on 
U.5U + l.ZU 


1 aa i 1 ^n 
1.00 + I.jU 


JUUU 


A^nn 
OjUU 


o &m~c no 


/JJ / 


1Z.J44 


noon 
U.Z8U 


V^A1 T \ir 

v joi Lyr 


1 OA i n S7 
1 .ZO + U.o / 


l.UZ + u. /z 


ozuu 


+ jUU 


Z. 1 / oJi-U / 


1 1 no i 

1 1UZ1 


J. loJ 


U.l JO 


Kl SCl 


1 AQ i A 


1 A7 i 1 no 
1.0/ + l.UZ 


7nnn 
/UUU 


a onn 
4SUU 


n chac no 


oqaa 
y3UU 


QA QAQ 


u.uyz3 


CL Aur 


1.35 + 2.24 


2.51 + 2.58 


6323 


9420 


i i mc r\H 
l.JUih-U/ 


10598 


45.080 


0.178 


u v^,cp 


1 Sfi 4- ^ 57 

1 .OU J .J / 


4.40 + 2.41 


407 5 


1 1215 


S 0Q7F-07 






n ^77 

\J.J I I 


USge 


1.99 + 5.45 


5.64 + 4.11 


5500 


12250 


2.036E-06 


20000 


49.185 


0.510 


SVCen 


8.56 + 6.05 


5.90 + 5.00 


14000 


23000 


1.626E-04 


37580 


211.238 


0.0297 


System 


M d + M g 


Rd + Rdisk 


^//.rf + T eff,g 


Tdisk f edge, disk 


dM/dt 


T 

± spot 


K 


D 




(M ) 


(Re) 


(K) 


(K) 


(Meyr- 1 ) 


(K) 






SW Cyg 


0.50 + 2.50 


4.30 + 3.44 


4891 + 9000 


6308 + 4968 


2.130E-07 


13060 


106.032 


0.587 


V356 Sgr 


3.00+ 11.00 


13.20 + 9.07 


8600 + 16500 


6174 + 4299 


4.442E-07 


17050 


659.326 


0.561 




4.25 + 14.1 


16.70+ 15.88 


13000 + 28000 


18279 + 8919 


3.440E-05 


22590 


135.355 


0.446 



Notes. The many digits in the numbers are only significant because they assure the internal consistency between the values of L acc , A, S and K. 
The first ten binaries are direct impact systems, the last three lines contain disk systems, with radii and temperatures of the gainer and its accretion 
disk defined separately. 



temperatures of both binary components. The error on the sum 
of the masses of a binary (horizontal axis of Fig. [3} is not fre- 
quently mentioned by the authors. In some other cases an error 
on the masses could be evaluated from various values found in 
the literature for a particular system. 

Moreover, the calculation of K uses the knowledge of the 
mass transfer rate (M* £0F ), the luminosity added to the gainer 
by the impinging mass (L fl ^), the surface of the hot spot (S spot ) 
and the hot spot temperature (T spot ). These observed quantities 
are theoretically interdependent. The quantity K can e.g. be cal- 
culated from Eq. {5} and/or from the relation K = (A Sy 1 . The 
quantity A can be determined from observation A = j^- and/or 
calculated from S and T spot . The quantity S can be directly ob- 
served and/or calculated from Eq. d3j. Erroneous evaluations of 
the observations can therefore be rejected when they violate the 
required internal consistency, as illustrated in Sect. A. 3 This is 
the main reason why only 13 binaries could be used in this anal- 
ysis. Errors on K (vertical axis of Fig. [3jl are obtained when the 
quantities mentioned above can be combined into models that 
respect every observed quantity within a reasonable range. The 
length of the vertical error bar then shows the largest and the 
smallest value that can be calculated for K within this frame- 
work. No error is shown in Fig. [3] when the error is too small or 
when not every quantity has been measured so that the internal 
consistency could not be completely explored. 



2.4. The distribution of mass ratios and orbital periods of 
Algols 

2.4.1. The observed distributions 



with an early B-type prim ary a t bi rth. When discussing initial 

M. 



conditions (e.g. in Sections 



2.2 



and 



3.1 



we use a = tt as def- 

1 Md 



inition of the mass ratio. However, since in Algols the gainer 
has become the most massive component, we use q-^-as the 
definition of the mass ratio of an Algol system. 



2.4.2. The calculated distributions 



Using the initial conditions of Sect. 2.2 and the data of our li 



brary with liberal binary evolutionary calculations yielding the 
durations that an Algol lives with well defined values of the or- 
bital period and mass ratio, Van Rensbergen et al. (120101 ) used a 
Monte Carlo simulation to calculate distributions of mass ratios 
and orbital periods of Algols that started RLOF during core hy- 
drogen burning of the donor (case A). This paper extends these 
calculations with Algols issuing from binaries that started RLOF 
during hydrogen shell burning of the donor (case B). We com- 
pare the results obtained under the assumption of strong tidal 
interaction with those obtained with weak tides. Although sig- 
nificant differences frequently occur between both assumptions, 
they yield a very similar overall distribution of mass ratios and 
orbital periods of Algols. 

Systems which become too faint with increasing distance to 
be included in the catalogs are removed by our Monte Carlo sim- 
ulation. Therefore we added the spatial distribution of binaries 
into the simulation so that they have either a spherically uniform 
distribution around the Sun or that they are uniformly distributed 
in a fiat structure such as Gould's belt. These two different as- 
sumptions yield also however no significant nor systematic dif- 
ferences. 



The distribution of orbital periods of Algols is well estab- 
lished. The distribution of mass ratios is however biased. Van 
Rensbergen et al. ([2008 ) used the mass ratio and orbital period 
distribution of Algols out of a sample of 303 observed systems. 
These systems were taken from the catalog of Budding et al. 
(2004) extended with semi-detached Algols from the catalog of 
Brancewicz et al. Q980). All these systems are issued from a bi- 
nary with a B-type primary at birth. A majority of 268 systems 
have a late B-type primary progenitor, leaving only 35 systems 



3. Results 

3.1. Calculated mass driven out of a close binary 

Using the stellar evolution code to calculate the evolution of bi- 
naries including those performing case B RLOF, we find that the 
amount of mass driven out of a close binary (AM) increases with 
increasing initial mass of the donor and with increasing initial 
orbital period. Binaries having an initial mass of the donor star 
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Fig. 3. Radiative efficiency of accretion K as a function of the total mass of a binary. Relation (|6| is shown by a thick line. The two 
thinner lines are at one standard deviation above and below the thick line. Errors and realistic estimates of inaccuracies are shown 
for individual cases. Filled squares indicate the positions of the ten direct impact systems whereas filled circles are used for the three 
systems with an accretion disk. 



below 6 M lose no or almost no matter during their evolution. 
This is because the mass transfer rate in these systems never ex- 
ceeds the critical value which is required in order for the combi- 
nation of spin-up and accretion luminosity to overcome the grav- 
itational potential of the system. Binaries with donors having ini- 
tial masses between 6 M Q and 9 M undergo a moderate amount 
of mass loss starting at larger initial orbital periods. Binaries with 
an initial mass of the donor star above 9 M suffer from severe 
mass loss during their evolution within this scenario. Although 
large individual differences may occur between binaries evolv- 
ing with strong, respectively weak, tidal interaction, there are no 
significant differences between the two groups. Figure [4] shows 
an example of the amount of mass lost by a binary with a 7 M Q 
primary at birth and evolving with strong tidal interaction. The 
grid of calculated cases covers initial mass ratios q = 0.4, 0.6 
and 0.9. At small initial orbital periods no mass is lost because 
RLOF is ignited before the donor has become large enough to 
transfer mass to the gainer at a rate larger than the critical one. 
At higher initial orbital periods most mass is lost in the q = 0.6 
case, sufficiently large to ensure a large enough separation (ac- 
cording to Kepler's third law) so that the mass transfer rate is 
higher than with q - 0.4. The mass lost in the case of q = 0.9 
is however smaller because such systems grow into near contact 
more easily, leading to a small geometrical factor D and thus to a 
fading accretion luminosity. The largest amounts of mass are lost 
in those cases where RLOF starts after exhaustion of hydrogen 
in the core of the donor star. This will be achieved in this case 
with initial orbital periods above ~ 4 days. The liberal calcula- 
tions in our catalog (Van Rensbergen et al. 2008 ) give also the 
amount of mass that is lost by many binaries with initial mass of 
the primary between 3 M Q and 15 M Q . 

3.2. Period distribution of Algols with a B-type primary at 
birth 

Figure [5] compares the observed orbital period distribution of 
Algols with a B-type primary at birth with theoretical simu- 
lations as obtained with the liberal assumptions. This distribu- 
tion is not affected by the different assumptions that were made 
about the tidal interaction. Algols observed with orbital periods 



of more than 15 days which were not reproduced sufficiently 
by considering the case A RLOF alone (Van Rensbergen et al. 
120101 ) are very well represented now that case B is included 
in the calculations. The maximum deviation in a one-sample 
Kolmogorov-Smirnov test decreases moreover from 0.115 ob- 
tained with the conservative simulation to the significantly better 
value of 0.067 when the liberal simulation is used. The observed 
period distribution of Algols is however only very closely re- 
produced by the liberal simulation for periods above 1 day. The 
number of ~ 9 % of Algols observed with orbital periods below 
1 day is largely underestimated by our simulation because these 
binaries tend to merge in our calculations. The fact that their 
number is very well observed could mean that binaries can sur- 
vive overcontact better than tolerated by our code (see discussion 
in Sect.|4|. 

3.3. The mass ratio distribution of Algols with a B-type 
primary at birth 

We used an observed distribution of Algols with a B type pri- 
mary at birth, having ~ 45 % in the high mass ratios (q e [0.4,1]). 
More specifically we find ~ 24 % with q e [0.4,0.55] and ~ 21 % 
wifhge [0.55,1]. The fact that the conservative binary evolution- 
ary code did not reproduce this observed mass ratio distribution 
for a sample of 303 Algols was the initial reason for us to start 
thinking about a liberal scenario. 

Van Rensbergen et al. ([2006 ) used /3 as a parameter taken 
equal to 0.5 during the entire binary evolution. They found that 
this a priori liberal assumption fits the observed mass ratio dis- 
tribution better than the conservative assumption does. In our 
subsequent papers we upgraded the quantity from a free constant 
parameter to a time dependent quantity /3(t) which is determined 
by the physical characteristics of the system. 

Conservative codes produce only ~ 12 % of Algols with a 
mass ratio q above 0.4. Including case B binary evolution to our 
statistics, this fraction rises up to ~ 17 %. Both results are shown 
in Fig. [6] It is clear that this limited enhancement is entirely due 
to those binaries with an initial primary mass above ~ 7 M Q be- 
cause in our scenario the less massive systems hardly lose any 
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Period (d) 



Fig. 4. Calculated amounts of mass lost by a binary with a 7M Q primary at birth undergoing strong tidal interaction during its 
evolution. The total amount of lost mass AM is on the vertical axis, the initial orbital period is on the horizontal axis. Calculated 
cases are indicated by filled circles. Binaries starting with an initial mass ratio q = 0.4 are connected by a long-dashed line, those 
with q = 0.6 by a dashed line and those with q = 0.9 by a solid line. 
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Fig. 5. Observed distribution (white) of orbital periods of 303 Algols as compared to liberal (black) simulation. 



mass during their evolution. Moreover only ~ 2 % of the Algols 
are reproduced with a mass ratio above 0.55, a number which is 
largely inferior to the ~ 21 % required by the observations. 



In Sect. 3. 1 we showed that binaries with an early B-type pri- 
mary at birth lose a large fraction of the transferred mass during 
their evolution. The short liberal era during which mass is lost 
(J3 < 1) occurs when the mass transfer rate is large, i.e. soon after 
the onset of RLOF when the binary is in its pre- or early- Algol 
stage. The mass ratio distribution as obtained with the liberal as- 
sumptions for binaries with an early B-type primary at birth will 
differ very much from the one obtained using the conservative 
binary evolutionary code. Hence we obtain ~ 33 % of Algols 
with q € [0.4,0.55] when we consider the liberal evolution of bi- 
naries with a mass of the initial primary above 7 M . Above 8 
M Q this fraction even rises up to ~ 45 %. The fraction with q e 
[0.55, 1 ] remains however always in the vicinity of 2 %, a number 
which is far below what is required by our set of observations. 



4. Discussion 

- Since more Algols are observed with orbital periods below 
one day than are produced by either evolutionary code, we 
propose that theoretical evolutionary calculations for bina- 
ries should allow for the survival of the system even after 
deep contact. This is similar to stating that such systems 
seem to survive a specific kind of common envelope evo- 
lution involving two main sequence stars. The nature of such 
evolution requires further exploration, as the physics are 
probably different from and more complex than the physics 
in the classical common envelope scheme as described by 
e.g. Webbink JT9841 and Nelemans & Tout (EOOBl 

- Although the theoretically calculated liberal mass ratio dis- 
tribution of Algols with an initial early B-type primary fits 
the observations better, the overall observed mass ratio dis- 
tribution of Algols still shows too many systems with high 
mass ratios. This could partly be due to theoretical mass 
transfer rates which are too large and too short living dur- 
ing the rapid era of mass transfer. We have compared mass 
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Fig. 6. Observed distribution of mass ratios of 303 Algols (white) as compared to conservative (grey) and liberal (black) simulation. 



transfer rates as obtained by our binary evolutionary code in 
the conservative mode with mass transfer rates for conserva- 
tive evolution as produced by previous authors and published 
by Kippenhahn & Weigert (11967b . Kippenhahn et al. (11967b 
and Paczyhski ( 1967a] [bj. The durations of our calculated 
phases of rapid mass transfer are very similar to those men- 
tioned above, whereas our peak values are somewhat lower 
because our stellar models are calculated with OPAL opac- 
ities (Rogers & Iglesias 119921 Iglesias et al. 119921 Iglesias 
& Rogers 1996), which were not available previously. Our 
calculated durations and peak values of mass transfer rates 
agree very well with those published by Nelson & Eggleton 
(2001) for a representative set of interacting binaries. The 
occurrence of many observed Algols with high mass ratios 
thus remains unexplained. 

Our scenario predicts only ~ 2 % of Algols with a mass ra- 
tio above 0.55. It is not certain that the discrepancy with ~ 
21 % of observed Algols with these large mass ratios is en- 
tirely due to fundamental failures of the theory. We would 
suggest that the determination of mass ratios from observa- 
tions is open for revision. Algols with a high mass ratio are 
present mainly when determined using main sequence char- 
acteristics for the gainer {qms)- Other methods generate an 
observed mass ratio distribution with much less high values 
(see Sect. 3 of Van Rensbergen et al.|2006). The validity of 
one method over another should be clearly decided before 
any conclusions can be attached to the difference between 
observed and predicted number of binaries in these high-q- 
bins. 

Call for observations: Much work has to be done in the fu- 
ture on the calibration of the efficiency by which mechanical 
energy is transformed into radiation as represented by rela- 
tion (J5J. This can only be done when data for Table[TJcan be 
refined and extended. The crucial quantity K is on one side 
determined by relation Q and on the other side by relation 

Relation Q uses the product of A > 1 with S « 1. The 
quantity A is determined from the knowledge of L acc and 
L ac M- Assuming that the geometry of the binary is per- 
fectly known, the value of L acc is determined with the (O- 
C) method (Sterken 2005) yielding the mass transfer rate 
only when the ephemeris of the eclipse changes linearly with 
time. The luminosity added to the gainer L a dd has been de- 



termined in only a few cases. The fractional surface has only 
been measured sometimes since small spot radii (< 10°) are 
hard to observe (Van Hamme et al. 12001b . Fortunately, we 
have at our disposal the theoretical relation as derived by 
Pringle (Eq. ([3])) which has no observational limitations and 
which is reliable within a factor of two for every case that 
could be checked in this paper. 

Relation <|3j uses quantities which are relatively well known 
such as the radius of the gainer and its effective tempera- 
ture. The uncertainty on the quantity L acc is discussed above. 
Measurements of the spot temperature T spot are crucial and 
are given in Sections 



A.l and A.2 



Binaries for which the value of K as determined from re- 
lation Q is not compatible within resonable limits on the 
uncertainties on S , L a dd, L acc , T spot with the value as deter- 
mined from relation <|3j have been removed from our sample. 
Although masses, radii and effective temperature of donor 
and gainer are usually well known, there is an explicit need 
for more and more reliable determinations of characteris- 
tics of the gainer star in semi-detached binaries. Observable 
quantities such as S , L a dd, L acc , T spot are needed to construct 
a better model which allows some binaries to lose mass into 
interstellar space. 

5. Conclusions 



- A catalog (Van Rensbergen et al.[2p08) with 561 conserva- 
tive and liberal evolutionary tracks is available at the Centre 
de Donnees Stellaires (CDS). Binaries with an initial pri- 
mary mass e [3,5] M Q are calculated in one mode only since 
they evolve conservatively. Binaries with an initial primary 
mass e [6,15] M Q are calculated in the liberal mode. Results 
for evolution with weak and strong tidal interaction are given 
separately. Conservative tracks are always added so that the 
reader is able to compare results of liberal (in two different 
tidal modes) and conservative evolution. 

- Mass impinging from the donor spins up the gainer and cre- 
ates a hot spot in its equatorial zone or at the edge of its ac- 
cretion disk. The combined energy of the enhanced rotation 
and increased radiation from the hot spot may exceed the 
binding energy of the system. We present the results of lib- 
eral evolutionary calculations for binaries with a B-type pri- 
mary at birth. A significant fraction of the transferred mass 
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is lost by the system in the case of binaries with an early B- 
type primary at birth. Systems with a late B-type primary at 
birth hardly lose any matter. 
- The observed distribution of orbital periods of Algols is well 
reproduced by our liberal theoretical calculations. Algols 
with an orbital period smaller than 1 day are however present 
for ~ 9 % in the observations whereas they are almost not 
created by our simulation. Our liberal theoretical calcula- 
tions reproduce the adopted observed distribution of mass 
ratios of Algol binaries less satisfactorily: Algols with mass 
ratios larger than 0.55 are underrepresented. Sect.|4]proposes 
possible additions to the present scenario to meet these dis- 
crepancies. 
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Appendix A: Data for 13 semi-detached binaries 

A. 1. Ten systems with direct impact 

Data on present masses and radii of the stars, mass transfer rates, 
temperature of the hot spot (T spot ), fractional size of the hot spot 
(S ), luminosity added by the hot spot (L„m) were taken for the 
ten systems with direct impact on the gainer's surface from the 
following references: 

- VW Cep: Masses, radii, effective temperatures of donor and 
gainer are from Pustylnik & Niarchos ( 20(30j. L a dd, T spot , 
R spot (and thus the quantity S) on the gainer are from the 
same authors. Values for P and thus of the mass transfer rate 
M are from Pribulla et al. (12000I ) and Devita et al. dl997l l. 

- AR Boo: Masses, radii, effective temperatures of donor and 
gainer and also values for T spot and M are from Lee et al. 
(I20Q91 

- CN And: Masses, radii, effective temperatures of donor and 
gainer and also values for T spot , S spo , and M are from Van 
Hamme et al. (2001) and references therein. A mean value 
of M was taken by considering also a somewhat different 
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value of M determined by Samec et al. 19981 1 and used by 
Qian (2002 ) in his analysis of near-contact binary systems. 
KZ Pav: Masses, radii, effective temperatures of donor and 
gainer and also values for L a dd and P are from Budding et al. 
(120011 1. We used the latter value of P in order to calculate M. 
V361 Lyr: Masses, radii, effective temperatures of donor and 
gainer are given by Hilditch (2001) and Yakut & Eggleton 
(2005 ). Hilditch (120011 ) also gives values for T spot , L a dd and 
M. This binary is thus very completely described. 
RT Scl: Masses, radii, effective temperatures of donor and 
gainer are given by Banks et al. (119901 ) and Pourbaix et 
al. ( 120041 ). Banks et al. (119901 ) also give a value for T spol . 
Values for M are given by Clausen & Gronbech d 19771 ) and 
Duerbeck & Karimie (119791 ). From this we calculated a mean 
value for M which is a factor of ~ 2 smaller than the value 
given by Rafert & Wilson 019841 ). 

CL Aur: All the data for this Algol binary, including a value 
for T spot on the gainer are given by Lee et al. ( 120 101 ). 
U Cep: Masses, radii, effective temperatures of donor and 
gainer are given by Budding et al. (2004). A value for 
P=+9.796 10~ 7 dy~ l was derived by ourselves from the 
"Eclipsing Binaries Minima Base" of Kundera, completed 
with the epochs of the most recent eclipses that can easily 
be found on the CDS website. From this we derived a value 
of M which is a factor of ~ 2 smaller than the evaluation by 
Pustylnik 019951 The value o f T spot e [20000, 40000] K was 
determined by Peters ( 120071 ) from the observations with the 
FUSE satellite. A refined value of T 



spot 



30000 K was then 
adjudged to the hot spot by the same author if one assumes 
that this HTAR (High Temperature Accretion Region) radi- 
ates as a black body. It is also interesting to notice that the 
gainer in this binary rotates 7.8 times faster than the synodic 
value (Dervisoglu et al. 120 101 ). 

U Sge: Masses, radii, effective temperatures of donor and 
gainer are given by Kempner & Richards 01999 ) and Vesper 
et al. (2001 ). Manzoori & Gozaliazl (2007 ) have determined 
a value for P and thus M and have seen a hot spot on the 
gainer's surface for which a range of temperatures has been 
given by Richards & Albright 019991 ). A value for L add has 
been determined by Albright & Richards 01 9961 ). This makes 
this binary also a well described system. 
SV Cen: Masses and radii of donor and gainer are given 
by Rucinski et al. (1992) because this author considers SV 
Cen not as a contact binary but as a semi-detached binary. 
The effective temperatures of both components are taken 
from the calibration for a B6.5III-donor and a B2V-gainer. 
Although Wilson & St arr 01976] ). Drechsel et al. 0T9g2"b and 
Herczeg & Drechsel 019851 ) consider SV Cen as an over- 
contact binary, many of their observations have been in- 
terpreted with SV Cen as a very active interactive binary. 
Herczeg & Drechsel ( 1985 ) report a large but variable value 
of P. Large mass transfer rates of ~ 4 10~ 4 M Q yr^ 1 and 
"a few times" 1O~ 4 M yr~' have been calculated by Wilson 
& Starr 019761 ) and Rucinski et al. 019921 ) respectively. 
Moreover a very large value for the hot spot temperature 
(Tspot ~ 10 5 K) and a very small value for the relative size 
of the hot spot (S < 0.01 ) have been reported by Drechsel et 
al. ([1982). In order to keep the values of M (and thus L acc ), 
A, S and K internally consistent and physically acceptable, 
we used in Table [ljfhe va lue of M = 1.608 10~ 4 M o yr 1 as 
given by Hilditch (2001 ), leading to a spot temperature of ~ 
37500 K which is well below the extremely large value of ~ 
10 5 K cited above. 



A.2. Three gainers with an accretion disk 

We have three disk systems in our sample. A hot spot temper- 
ature is available for j3 Lyr only. The quantity L a dd which is 
needed to calculate the hot spot temperature is calculated with: 



Lgdd 



spot 



/ 7"* _ T"4 "\ 

spot 

5770 4 



(A.l) 



The hot spot is now on the edge of the disk and its surface 
has never been measured. Expression ([3]) was derived for the cal- 
culation of the size of a hot spot on the surface of the gainer and 
can thus hardly be used in this case. Using this expression never- 
theless in the case of f3 Lyr we find an amount L acc that exceeds 
the calculated accretion luminosity by a factor of three. In or- 
der to obtain the added luminosity we have therefore divided the 
surface area of the impact zone as given by Pringle (1985]) by 
the same factor of three for the three cases in this subsection. In 
equation (A.l i we used the effective temperature of the gainer 



for the two systems where the hot spot is located nearby the sur- 
face of the gainer (SW Cyg and V356 Sgr) and we used the tem- 
perature in the middle of the accretion disk where the spot is 
visible as a source of bipolar jets for /3 Lyr (Linnell 120021 ). In 
order to calculate K we replaced in relation (pMC by Rd ls k and 



T e ff,g by T e d g e,disk, using relations (A.2) and (A.3 i 



The characteristic disk temperature and the temperature dis- 
tribution from the stellar surface to the edge of the disk were cal- 
culated from the model of Bath & Pringle 01981I ). Home ([1985 ) 
developed a method of eclipse mapping of accretion disks with 
this model and Rutten et al. ( 1992 ) applied the latter method to 
construct models of accretion disks for six cataclysmic variables. 
The values of Tdisk an d T d isk.edge mentioned in Table [T] were thus 
calculated with: 



'disk ~ 
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= 478074 K 
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(A.2) 



with ctr = 5.67 10 5 erg cm 2 K 4 s 1 the Stefan-Boltzmann 
constant and 



R g 

T(r) = Tdisk I — 



(A.3) 



For the three systems with an accretion disk around the 
gainer, data on present masses and radii of the stars, mass trans- 
fer rates, temperatures of the hot spot (T spot ) and added lumi- 
nosities through the hot spot (L a dd) were taken from the follow- 
ing references: 

- SW Cyg: Masses, radii and effective temperatures can be 
found in the catalog of Budding (120041 ) and Richards & 
Albright 019991 ). Values for P (and thus M) are from Qian 
et al. 020021 ). Albright & Richards 019961 ) argue that the lu- 
minosity added to the similar system TT Hya by accretion is 
five times as large as compared to SW Cyg. Since Peters & 
Polidan (2004) establish a hot spot temperature of 17000 K 
in TT Hya (more than 60 % hotter than the underlying effec- 
tive temperature of the gainer), we conclude that the obtained 
hot spot temperature of ~ 13000 K for SW Cyg is reasonable. 

- V356 Sgr: Masses and radii of donor and gainer were found 
by Peters & Polidan ( 1998). These authors also mention the 
existence of very hot circumstellar material around this bi- 
nary and of a hot spot concentrated near the photosphere of 
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the B-type gainer. Effective temperatures of donor and gainer 
are from Polidan ( 1989 ). We derived a value of M from the 
quadratic term in the ephemeris of eclipses as given by Hall 
et al. ( 1981 ). This value is in good agreement with the mass 
transfer rate as given by Polidan ( 1989 ). It is pointed out by 
Simon ( 1999 ) and Dervisoglu et al. (120101 1 that the gainer ro- 
tates at the the large equatorial velocity of ~ 212 km s . This 
is however far below the critical velocity of ~ 612 km s~ l 
so that in view of the Roche model (e.g. Maeder 2005) the 
gainer remains practically spherical, and the observed UV- 
excess cannot be due to the oblateness of the gainer with 
polar regions much brighter than the equatorial ones, as sug- 
gested by Polidan (f!M9"l) . 

- p Lyr: Masses and radii of donor and gainer are from Simon 
(1999). Effective temperatures are from Harmanec (2002). 
We derived a mass transfer rate from the values of P given 
by Simon ( [19991 and Ak and al. (120071 . Harmanec (120021 
reveals a minimum temperature for a hot spot of 20000 K. 
Linnell (2002) proposes that the kinetic energy of the mass 
transfer stream is converted into heat and radiation at a radius 
equal to half of the radius of the accretion disk. Schmidtt et 
al. ([2009 ) observed strong Ha originating from the disk, pos- 
sibly at the point of impact of the gas in the disk. The disk 
temperature of 8500 K as mentioned by Harmanec (2002) 
is very similar to the temperature we find at the edge of the 
opaque accretion disk of ji Lyr (~ 8920 K). The masses of 
the components as derived by Zhao et al. (120081) are some- 
what lower than those mentioned by Simon (1999). They 
yield somewhat lower values of the disk temperatures and 
somewhat larger values for K. They are included in the error 
bars in Fig. [3] 

A3. Systems removed from the sample 

We had to exclude many binaries from our statistical analysis. 
Some systems are indeed extensively studied with one crucial 
quantity that was however never determined from the observa- 
tions. For other systems we found observational values of phys- 
ical quantities which were contradictory. We mention a few ex- 
amples of systems which were for the time being excluded from 
our statistics. 

- TT Hya: Many data on this Algol system with an accre- 
tion disk are known, including the value of 17000 K for the 
hot spot temperature by Peters & Polidan (2004). The mass 
transfer rate of > 2 10~ 10 M Q yr~ l as found by Miller et al. 
(2007 ) is however five orders of magnitude too low in order 
to obtain such a high hot spot temperature. The quoted value 
of the mass transfer rate would yield a hot spot with a tem- 
perature of 1000 K at most: only a little bit higher than the 
temperature at the edge of the accretion disk. 

- j6 Per: The same can be said about the one and only real 
Algol, /3 Per. Images of accretion structures containing a 
stable emission source on the gainer were obtained with 
Doppler tomography by Richards (2004). A high tempera- 
ture region (T « 10 4 - 10 5 K) surrounding the gainer has been 
derived from the observation of high temperature resonance 
emission lines by Cugier & Molaro ( 119841 ) and Brandi et al. 
(1989 ). Blondin et al. d 19951 ) found that the gas stream to- 
wards the gainer is heated up to a temperature of 2 to 4 10 6 K 
for this direct impact system. Such high temperatures can 
not be produced by a mass transfer rate as low as 0.4 to 2 
1O~ U M yr' 1 as proposed by Pustylnik ( 119951) . A hot spot 
with a temperature of ~ 10 5 K can only be created through a 



mass transfer rate which is seven orders of magnitude larger. 
A mass transfer rate as low as ~ 1O~ U M yr~ l would raise 
the photospheric temperature of the gainer with less than 1 K 
according to the scenario used to construct Table [T] 
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